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1. Introduction 

The relevance of weak structural disorder for the critical behaviour of pure magnetic 
systems has been in the focus of attention of researchers since a long time Depending 
on the structure of a magnet disorder can enter as random fields, random anisotropics, 
random bond or site dilution. Our special interest here is in diluted systems. Concerning 
the statics of such systems one of the main insights comes from the Harris criterion |2] 
stating that the static asymptotic critical behaviour remains unchanged, if the specific 
heat of the pure system does not diverge P^. Since the borderline value ric between a 
diverging and nondiverging specific heat in space dimension c? = 3 lies between order 
parameter (OP) dimension n = 1 (Ising model) and n = 2 (XY model) |3] only diluted 
Ising systems show a new critical behaviour, which is characterized by a non-diverging 
specific heat p. 
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Systems with the same static critical behaviour may have different dynamic critical 
properties. In the vicinity of a critical point slow processes play an important role, 
therefore in critical dynamics the main objects are the slow modes. Apart from the 
order parameter these are the conserved quantities, which coupled to the OP. The 
simplest dynamical model - model C |H] - includes these couplings only statically. For 
such a case the Harris criterion also has consequences. It leads to the conclusion that 
the coupling of conserved quantities to the OP in disordered systems is of no relevance 
[01 13 IE] • The argumentation of this statement is based on the fact [3 El HD] that 
coupling to a conserved density for the relaxational model is relevant only if the specific 
heat diverges. Diluted models always have nondiverging specific heat and therefore their 
coupling beetwen a conserved density and OP is irrelevant. As a consequence most of 
the studies of critical dynamics considered only relaxational dynamics of diluted Ising 
systems HH [121 UHl d • 

However the conclusion concerning the irrelevance of the conserved density is 
made only for the asymptotic properties of the model. Meanwhile it became 
clear that in diluted magnets most of the experiments and simulations are in the 
nonasymptotic region jH [13 [lEl [HI [IHl [IHI • Within the standard tool of investigation 
of critical phenomena, namely the field-theoretical renormalization group (RG), the 
nonasymptotic critical behaviour is described by the flow of the model parameters (static 
and dynamic) leading to effective critical behaviour rather than to the asymptotic one 
characterized by power laws with universal exponents determined by the fixed point 
properties of the parameters. 

Thus we consider the dynamics of a diluted magnetic system where the OP relaxes 
with a relaxation rate F and the second conserved density shows diffusive behaviour 
with a diffusion rate A. The ratio of the time scales = F/A is the dynamic parameter 
of model C and takes on different fixed point values depending on the dimension n of 
the OP and spatial dimension d. Only recently the correct field theoretic RG functions 
of model C have been calculated Ej. It turned out that model C itself has a very 
slow transient at least for the cases n = 1, 2 in two-loop order. 

Another artifact present in the RG analysis of diluted magnets within one loop 
order is the degeneracy of the static /^-functions for n = 1. As it is well established 
now it leads to a ^/A — d rather than a 4 — (i-expansion [2111 [^ • This makes a two- 
loop treatment inevitable. Moreover the borderline where the specific heat exponent a 
changes in two-loop order is shifted from = 4 to below n = 2. 

Thus two-loop calculations are necessary (i) to specify the shift of the stability 
regions and (ii) to calculate quantitative values of the effective exponents. 

All the above arguments serve as a reason to consider the disordered model C 
critical dynamics by the field-theoretical RG approach within two-loop approximation. 
Being more elaborated technically this approximation should lead to reliable numerical 
results which as previous experience shows should not be changed essentially by further 
increase of the perturbation theory order. Some of our results are briefly summarized 
in Ref. [SUj, and here we give a more thorough derivation of the dynamic RG function 
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and analyze the diluted model C effective critical dynamics. 

The set-up of the paper is the following: in the next Section we define the dynamical 
model for disordered magnets, its renormalization and general relations for the field 
theoretic functions. Results in one-loop order are collected in Section El Two-loop 
results are presented in Section |3] Conclusions are given in Section IHl 



2. Model and renormalization 



2.1. Model equations 

The object of our analysis consists in a dynamical model for quenched disordered 
magnets, namely model C in the classification of Ref.j22|- Model C describes the 
dynamics of a nonconserved OP (p^ which is coupled to a conserved density uiq (in most 
cases the energy density). The secondary density has to be taken into account since it 
is also a slow density showing critical slowing down near the phase transition. The OP 
is assumed to be an n-component vector, while the density mo is a scalar quantity. The 
structure of the equations |221 IHI of motion is not changed by the presence of disorder. 
They read: 

^— =-r- + e^^ t = l...n, (1) 

dt dipifi 

— = AmV^TT \- 9m- (2) 

ot oniQ 

The OP relaxes to equilibrium with the relaxation rate P and the conserved density 
mo diffuses with the diffusion rate Am- The stochastic forces in ((H), Q satisfy the 
Einstein relations: 

< e^X^,t)e^^{x\t') > = 2t6ix - x')6it - t')5i,, (3) 

< em{x, t)em{x', t')> = - 2AmV'5(x - x')5{t - t')5i,. (4) 

In order to define model C in the presence of structural disorder we write the static 
functional Ti describing the behaviour of a disordered magnetic system in equilibrium: 

n= f d'^x{^rVo|' + nx)|^or + ^X^(V¥..,o)' + 

" l<^o|^ + ^ctmW^o^ + 7;^mmo\0o\'^ - hmmo], (5) 



where V{x) is an impurity potential which introduces disorder to the system, and d 
is the spatial dimension. The functional (0) contains a coupling 7^ to the secondary 
density which can be integrated out. Thus static critical properties described by the 
functional © are equivalent to those of a Ginzburg-Landau- Wilson (GLW) functional: 



n 



I rf'^x{ifi^or+v^(x)i^or+^x^(v^.o)^+|i^or}. (6) 
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The parameters r and u are related to f, m, a^, 7^ and hm by 

r = r H , u = u — 3 — , (7) 

r is proportional to the temperature distance from the mean field critical temperature, 
u is positive. 

The properties of the random potential V{x) are governed by a Gaussian 
distribution 



/ d'^xV{ 



X 



Vv = Mv exp ^-^ J ■ 

with positive width A, which is proportional to the concentration of non- magnetic 
impurities. 

For investigation of statics in the further procedure one has to take into account 
that the disorder is quenched. Averaging the free energy of a disordered system over 
the distribution (jH)) with application of the replica trick [221 one ends up with an 
effective static functional containing new terms determined by disorder j2Tj. Then 
critical properties are studied on the basis of long-distance properties of the effective 
functional. 

The procedure for the critical dynamics may be different. We will treat the critical 
dynamics of the disordered models within the field theoretical RG method based on 
the Bausch-Janssen- Wagner approach |21I, where the appropriate Lagrangians of the 
models are studied. Therefore we have to obtain the Lagrangian for our model on the 
basis of the model equations (0)-®- Then after averaging over the random potential 
(jHl) we get new terms determined by the disorder in the Lagrangian. In this case it is 
not necessary to apply the replica method Results can be written in the form (for 
details see [Appendix A| ) C = + Lmt with the Gaussian (unperturbed) part: 



d'^xdtl -f ^ (^o,i<^o,*+ 5Z <^o,* ( (r-V' 
i=i j=i ^ 

+A„moV^mo + mo - ctmAmV^^ moj, (9) 
and an interaction part: 

l^int= I ci'^xrft^ l^fwc^o.iV'o.iX^'^O'i'^oj- I dt'^2t'^A(fo^i{x,t)^o^i{x,t) 

X(fo,j{x,t')Lpoj{x,t') + f7mmo<^o,jV5o,i-^Am7m"^oVVo,jV5o,i| , (10) 

with new auxiliary response fields <^o,i- The ratio u/ A determines the degree of disorder 
in the system. 

Investigating the long-distance and long-time properties of the theory with the 
effective Lagrangian C we apply Feynman diagram techniques in order to get dynamical 
vertex functions. Details of the calculation of the dynamic vertex function F^,^ 



are given in Appendix A The renormalization of vertex functions leads to the RG 
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functions, describing the critical dynamics of our model. We use the minimal subtraction 
scheme with dimensional regularization to calculate these functions. General relations 
are considered in the next subsections. The details concerning the renormalization 
procedure are given in [Appendix B 



2.2. RG functions 

From the renormalizing factors introduced in Appendix B we define the ^-functions 
which describe the critical behaviour of our model 

(") 

here {a.i\ = {n, A,7,r,A} represents the set of renormalized parameters for the 
disordered model C and /i is the scale. The notation denotes any density {ip, (p, 
m, m) or any model parameter ctj. 

A special case is the additive renormalization A^2 of the specific heat in the GLW- 
model. It leads to an additional RG function 

BAu, A) = fi'Zl,fi4- (Zjf^-'A^.) , (12) 



which appears in the static and dynamic ^-functions. 

Since in statics the secondary field can be integrated out one finds relations to the 
static Z-factors appearing in model A (see Appendix B). The relations ()B.9|) and ()B.7|1 



lead then to 

C^(m,A,7) = 2C™(m,A,7) + C^(m,A) + C^2(m,A), (13) 

and 

Cm{u,A,-f) = ^-f^B^2{u,A). (14) 

Eliminating Cm by inserting Eq. (fT^ into Eq. (fT^ one gets: 

C^(m,A,7) = 725^2(n,A) + g(?i,A) + g2(M,A). (15) 

The (^"functions for the kinetic coefficients F and A, which are obtained by 
substituting relations ()B.5|) and ()B.13|) into the definition (fTT|) . read 

Cr({n, A, 7, w})=-^C^{{u, A, 7, w})+^CAu, A), (16) 

CA(n,7)=7'5^2(n,A). (17) 

The last relation is obtained taking into account Eq. ()14|) . 

In order to investigate the dynamical fixed points of model C it is convenient to 
introduce the time scale ratio w = r/X. The corresponding (^-function for this dynamical 
parameter using the relations (fTH|) and (fT7|) reads: 

C^(m, a, 7, w) = ^(^{u, A) - ^C>^(m, a, 7, w) - 7^S^2 (m, A) . (18) 

The static critical properties of our model are described by the fiow equations: 

l^ = P^{u,A), l^ = p^(u,A), /^=/3,(n,A,7), (19) 
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with the flow parameter ^ and /5-functions generally defined as 

= a^[-C^ + Ca.({«i})], (20) 

where q are the engineering dimensions of the corresponding parameters aj. Following 
Eq. ()2Up we can rewrite the static /^-functions: 

(3u{u, A) = n {-e - 2C^{u, A) + („(«, A)) , (21) 
/3a (w, A) = A {-e - 2C^(m, A) + Ca(w, A)) , (22) 

(3^{u, A, 7) = 7 - I - g(n, A) - C,n(n, A, 7) + C^(u, A, 7)^ . (23) 

Note that the fiow equations for u and A are independent from those for the other model 
parameters. They are the same as for the diluted n-vector model. 

The dynamical critical properties are determined by the fiow equation for the time 
scale ratio which reads 

=/5«,(m,A,7,^i;), (24) 

where the /3-function for w, according to Eq. (plj) . is defined as 

/?^(m, A,7,w) = ^^(m, A,7,w), (25) 

since F, A and their ratio w are dimensionless parameters. 

Using the relations between ^-functions for (3^ and I3w one obtains: 



(3^{u, A, 7) = 7 ( -- + CA^^ ^) + o^'^v^("' ^) 



l[-^^ + \l'BAu.^)]. (26) 



(3^{u, A, ^,w) = w Qc^(m, A) - ^C^(m, a, 7, w) - -f'BAu, A) j . (27) 

where A = e — 2(^2(u, A) equals the ratio of coupling dependent functions of the critical 
exponents of the heat capacity a{u, A) and the correlation length z/(m, A). 

2.3. Asymptotic properties 

The common zeroes of the /3-functions (|7 H) . (|^ . (j^ . (j?rj) define the fixed point (FP) 
values: {a*} = {u*, A*, 7*, w*}. The zeroes of the functions Pu and /5a can be obtained 
independent from the other /^-functions. For each pair of FPs {u*, A*} one obtains two 
values of 7* from /3^: 

7 =0 and 7 = = 28 

B^2{u*,A*) uB^2{u*,A*) 

where a and u are static heat capacity and correlation length critical exponent calculated 
at the corresponding FP {u*, A*}. Then the results for the static FPs are inserted into 
Piu in order to find the corresponding values of w*. 
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The relevant FP corresponding to a critical point of the system has to be (i) 
accessible from the physical initial conditions and (ii) stable. The stability of a FP is 
defined by the eigenvalues uJi of the matrix dj3j/dai. If the real parts of all uoi calculated 
at some FP {a*} are positive then the FP {a*} is stable and the flow of the system of 
differential equations (fTn|l and ()24j) is attracted to this FP in the limit £ — ^> 0. 

From the structure of the stability matrix we conclude, that the stability of any FP 
with respect to the parameters 7 and w is determined solely by the derivatives of the 
corresponding /3-f unctions: 

Moreover using ()26|) we can write: 



^. = -^^^%^ + ^7X^KA), (30) 
which at the FP a* leads to: 

^7U..: = for 7*^ = 0, (31) 

^7L.=„* =^ for 7*V0. (32) 

Therefore stability with respect to the parameter 7 is determined by the sign of a. For 
a system with non-diverging heat capacity (a < 0) at the critical point, 7* = is the 
stable FP. Diluted magnets we consider here always have a < p. This leads to the 
conclusion that in the asymptotic region the secondary density decouples from the OP 

PHIIH]. 

The critical exponents are defined by the FP values of the ^-functions. For instance, 
the asymptotic dynamical critical exponent z (at the stable FP) is expressed in the 
following way: 

^ = 2 + Cr(w*,A*,7*,«;*), (33) 

while its effective counterpart in the non-asymptotic region is defined by the solution of 
flow equations (|TT?|l and (j211) as 

Zefr = 2 + Cr(M(£),A(£),7(£),w(£)). (34) 

In the limit £ the effective exponents reach their asymptotic values. 



3. Results in one-loop order 

Although the one-loop order results have the drawbacks mentioned in the introduction 
one gets a quahtatively correct picture of the effects of disorder on the dynamics. 
Therefore we discuss this case in more detail. 

We are interested in the dynamical properties of disordered model C within the 
first nontrivial order of expansion in renormalized couplings, that is the one-loop 
order. The static functions (3u and j3/\ for disordered systems are known in higher 
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Table 1. One-loop FPs of model C with disorder and their stability for e = 4 — d > 0. 
Note that 7*^ has to be positive for the existence of a real FP. . 



FP 


u* 


A* 




P* 


stability 


G 











0< p* <l 


unstable 















G' 








n/2 


unstable 










1 




P 


6£ 

n+8 


n 
u 





< p* < 1 


marginal for n > 4 








2 (4— n)£ 
n n+8 







P' 






n/2 


unstable 










1 




U 





3£ 

4 








unstable 










1 















unstable 


U' 






€_ 

n 


3n/4 


stable for n < 4/3 










1 


stable for n > 4/3 


M 


3£ 


(4-n)£ 








stable for 1 < ?^ < 4 


2(ri-l) 


32(ri-l) 










1 


unstable 













unstable 


M' 






2 n-4 ^ 
n 4(n-l)'' 


3n/4 


ustable 










1 


unstable 



loop approximation jzoj. However taking the same order for statics as for dynamics, we 

restrict the expressions for j3u and /5a to the one-loop approximation: 

/ n + 8 \ 
Piu = ui-e + —^u - 24A j , (35) 

= a - 16A + ^^w) • (36) 

Note that the region of physical relevance of the couplings u, A for diluted magnets 
is restricted by m > 0, A > 0. The other is obtained using the static one-loop 
C-f unction: 

A, = 7(-2-4AH-— .+ -2^j, (37) 

For obtaining the dynamic function (5^^ we should first calculate the renormalizing factor 
Z(p. From the one-loop part of Zq, (see [Appendix A ) following formula (fTT|) we can obtain 
the C(^-function. Using obtained result together with one-loop static ^-functions we get 
in the following form: 

,._(4.,,=^-2^). (38) 

Setting the right hand side of Eqs. (jH^jl - ljHHj) equal to zero we obtain the system of 
equations for the FPs. The structure of the functions (5u and /5a leads to the existence of 
four FPs: the Gaussian FP G {m* = 0, A* = 0} , the 'polymer' FP U {m* = 0, A* ^ 0}, 
the FP of the pure system P {u* 7^ 0, A* = 0} and the mixed FP M {v* 7^ 0, A* 7^ 0}. 
Depending on the FP values u* and A* values for the genuine model C parameters can 
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be found. As it was already mentioned two values 7*^ = or 7*^ 7^ correspond to 
each FP, this doubles the number of fixed points, see Table ^ We note that the FPs 
with nonzero 7*^ are indicated by a prime '. We do not consider further the FPs U and 
U' since these FPs lie outside the physical region of positive values of A. 

Among the rest of the FPs only one is stable depending on the OP dimension n. 
For n > 4 FP P is stable while for 1 < < 4 FP M is stable. Thus the one-loop value 
of marginal dimension nc{e) = 4 defines the borderline where a = 0. In higher- loop 
order the FPs picture is not changed apart from a change of the borderline function 
nc{e). Estimates obtained on the basis of six-loop order calculations 3] give definitely 
ric < 2 at d = 3. 

In order to discuss the dynamical FPs it turns out to be useful to introduce the 
parameter p = w/{l + w) which maps w and its FPs on a finite region of the parameter 
space. Then instead of the flow equation the flow equation for p arises: 



Setting the right side of ()40|) to zero and using the FP values from Table ^ the 
values of p* can be found. They are shown in Table C] as well. 

Each non-zero solution of 7*^ leads according to ()39|) to three dynamical FPs: either 
with p* = (i.e. w*=0), p*=l (i.e. w*=oo) or < p* < 1 correspondingly. For the FP 
which corresponds to 7*^ = the situation is the following. For the FP M only two 
dynamical FPs with p* = and p* = 1 exist, while for the FPs P and G any value of p* 
between and 1 is allowed. Checking the stability of these FPs we see that for 1 < < 4 
only FP M with p* = is stable. The corresponding one-loop asymptotic value of the 
dynamical critical exponent in this case coincides with the one-loop result for the pure 
relaxational model (model A) with disorder: z = 2 + g'^~"|^ |1H |7j . For n > 4 FP P 
is marginal, all other FPs are unstable. Formally flow starting from the initial values 
with Pi ends up at some FP P with p*(pi) depending on the initial value p,. Thus one 
gets a whole line of FPs at u* 7^ 0, A* = 7* = 0. In this case the one-loop result for 
z coincides with the conventional theory value z = 2. Since in one-loop approximation 
nc{e) = 4 the FPs M and P determine the critical behaviour of the disordered model 
C for n < 4 and n > 4 respectively. 

4. Two-loop results 

Static RG functions are already known in high-order approximations within different 
renormalization schemes (for references see e. g. Ref jT]). Two-loop expressions for 
functions (3^, f3^ within the minimal subtraction scheme can be obtained in the 
replica limit from results of Ref. |3I] and they read: 



where according to 

Pp{u, 7, P) = p{p - l)Ct«(M, 7, p) ■ 



(39) 



(40) 




(41) 
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m + 2 



lA 



m + 2 r, 



= + 7 



36 
Sfm + 2 



3 

5(m + 2) 



-mA - 168A^ 



uA - 20A^ 



10 

(42) 
(43) 



72 3 

The dynamical function (3p is given by formulas ()18|1 and ()4()j) . where functions 
and are needed also for the calculation of the dynamical critical exponent is 
the same as in one-loop approximation. can be obtained in the replica limit from 
results of Ref . [SI] : 



m 



m 



-uA - 4A^ 



(44) 



72 3 

For the calculation of (^p we use the RG scheme described in detail in Appendices A and 
B. We obtain from the two-loop value of ()B.14|) , according to the definition (fTT|) 
it reads: 

'4 1 \ ^ w 
12 



G = -8A 



n + 2 , 
3 uA 



,0 n + 2 ^ 
44A2 



6 



In 



1+w 



+ 



n + 2 

w In 
T+w V2 

w 



M 1 - 3 In ■ 



w 



3(n+2) 



-I2A7 



l+w 



1+w 



+ 4A72 



3 



l+2w, il+w) 



In 



w 



l+w 



w In 



w 



l+w l+2w 
-31n(l + u;) - 



w 



7 



w 



l+w 
In w 



(45) 



l+w l+w 
Setting the coupling A equal to zero in (jl5|) . one recovers the two-loop result for pure 
model C |2H] , while expression with 7 = corresponds to the function C,g, of model 
A with dilution, which was extensively studied IT^ IT^ ITl] . The 7^A-term is an 
intrinsic contribution of the disordered model C. 

Two complementary ways are known for the analysis of the FP equations. The 
first one is the ^-expansion j^, while the second one consists in fixing e and solving 
equations for the FP numerically |SS1- As it is known from statics for diluted systems, 
expansions in e do not give reliable numerical estimates for critical exponents. For 
instance for disordered Ising magnets the e-expansion technique leads in fact to a \/e- 
expansion [201 12I|; that does not give trustable results for e = 1 |2E1- Therefore we 
follow the second way of analysis working directly in the d = ?> dimensional space. 

The series for RG functions are known to be asymptotic at best. For instance, 
no FP of the fourth-order couplings is obtained without application of resummation 
for diluted systems in the two-loop approximation jl]. Therefore for static RG 
functions it is standard to apply resummation technique in order to obtain reliable 
results. We use here the Pade-Borel resummation procedure |Mj for the functions 
/5a/A,/5^/7- 7 V4}. 



4.1. Fixed points and their stability 



We can analyze the RG functions (3u and /?a independently from the other ones. The 
results for these static functions are already known pp. The outcome of the two loop 
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approximation qualitatively repeats the results of the one-loop approximations, namely, 
the Gaussian FP G, the FP of the pure system P and the mixed FP M are found. 
However, contrary to the one-loop results, the FP M is defined also for the Ising case 
n = 1. 

Among the remaining FPs only one is stable depending on the order parameter 
dimension n. For n < n^, the FP M is the stable one, at n = FPs M and P change 
stability, and for n > FP P becomes stable. For the stability boundary at ci = 3 we 
find Tic = 1-55 which is somewhat smaller than higher- loop order estimations jl]. 

From the structure of two-loop (PHj) one concludes that all FPs described above 
exist for 7* = 0. Analyzing jS^ we find that only the Gaussian FP G and the FP P 
have counterparts with positive non-zero 7*: FP G' for all n and FP P' for n < 
respectively. The FP values of the model parameters and the dynamical exponents for 
all FPs are listed in Table El for the numbers of OP components n = 1, 2, 3. 

At FP G a line of FP for the time scale ratio < p < 1 is obtained for any n. 
For the FPs with 7* = only two values p* = or p* = 1 are found, while for the 
FPs P' and G' one obtains three solutions p* = 0, p* = 1 and a non-zero solution with 
< p* < 1. All FPs are given in Table El We denote the FPs with p* = 1 by subscript 
1, whereas FPs with p* = have no subscript. Since FP P' with non-zero solution 
< p* < 1 corresponds to the FP of pure model C, we denote it by C, and the FP G' 
with non-zero p* we denote by G',^. 

Analyzing the stability exponents u)^ and u)p we see that only FP M is stable for 
n = 1 and FP P is stable for n = 2, 3. That means that the model A critical behaviour 
is reached in the asymptotics in any case (diluted model A universality class for n < ric 
and pure model A universality class for n > ric). As a consequence the secondary density 
is for all n asymptotically decoupled and formally the dynamical critical exponent takes 
its Van Hove value 2;^ = 2, as expected due to arguments of Ref. [HllZl- 

The stability exponents are listed in Table El They determine the "velocity" of the 
RG flows in the FPs vicinity. A small value of a stability exponent calculated at a stable 
FP means slow approach to this FP in the corresponding direction. For example, for 
FP M at = 1 one has uj^ = 0.002, which leads to slow approach of 7 to its FP value 
7* = 0. 

In Table |21 we give the numerical values of the asymptotic dynamical critical 
exponents z calculated for all FPs. If the flow from the initial values of the couplings 
passes near one of these FPs one may observe an effective critical behaviour governed 
by the values of the critical exponents corresponding to that FP. 

4-2. Flows and effective exponents 

The nonasymptotic behaviour is described by the flow of the static couplings and the 
dynamic parameter under renormalization. It can be obtained solving the flow equations 
(fTTI|l and (jnni)- Using the solution we get the behaviour ofthe effective critical exponent 
Zc?i with continuous change of the flow parameter. 
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Table 2. Two-loop values for the FPs of disordered model C. 



n 


FP 


u* 


A* 


7* 


P* 


z 


Vn 


G 











< p* < 1 


2 


n = 1 


G' 

G'c 
G'l 
P 

Pi 
P' 
C 

Pi 
M 
Ml 






1.3146 
1.3146 
1.3146 
1.3146 
1.3146 
1.6330 
1.6330 











0.0209 
0.0209 


1.4142 
1.4142 
1.4142 





0.4582 
0.4582 
0.4582 








0.3446 
1 

1 


0.2664 
1 

1 


2 
3 

oo 
2.052 
2.052 
2.052 
2.105 
2.052 
2.139 
2.139 




G' 

G'c 
Gi 
P 
Pi 






1.1415 
1.1415 









1 
1 
1 






0.6106 
1 

1 


2 
3 

oo 
2.053 
2.053 


71 = 3 


G' 

G'c 
Gi 
P 

Pi 






1.0016 
1.0016 









0.8165 
0.8165 
0.8165 








0.7993 
1 

1 


2 
3 

oo 
2.051 
2.051 



Table 3. Stability exponents for the FPs given in Table [21 



n 


FP 




LUA 


LU-y 


UJp 


Vn 


G 


-1 


-1 


-0.5 





n = I 


G' 


-1 


-1 


1 


-1 




G'c 


-1 


-1 


1 


0.976 




Gi 


-1 


-1 


1 


-oo 




P 


0.566 


-0.105 


-0.053 


0.052 




Pi 


0.566 


-0.105 


-0.053 


-0.052 




P' 


0.566 


-0.105 


0.105 


-0.053 




C 


0.566 


-0.105 


0.105 


0.041 




Pi 


0.566 


-0.105 


0.105 


-oo 




M 


0.494 


0.194 


0.002 


0.139 




Ml 


0.494 


0.194 


0.002 


-0.139 


n = 2 


G' 


-1 


-1 


1 


-1 




G'c 


-1 


-1 


1 


0.745 




Gi 


-1 


-1 


1 


-oo 




P 


0.581 


0.078 


0.039 


0.053 




Pi 


0.582 


0.078 


0.039 


-0.053 


n = 3 


G' 


-1 


-1 


1 


-1 




G'c 


-1 


-1 


1 


0.522 




Gi 


-1 


-1 


1 


-oo 




P 


0.597 


0.222 


0.111 


0.051 




Pi 


0.597 


0.222 


0.111 


-0.051 
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First we consider the case of Ising spins n = 1. We present static flows in the 
subspace m — A — 7 as well as a "projection" of the dynamic flows on the subspace of 
the model parameters m — A — p in Fig. ^ All flows are obtained for initial values with 
7(^0) = 0.1. Flow equation solutions presented in Fig. [T]are obtained for different ratios 
A(£o)/m(^o) as well as different values of p(^o)- Flows starting from initial conditions 
with small ratio A(£o)/w(^o) (curves a,c,e) correspond to a small degree of disorder. 
They are influenced by the unstable FPs P, P', C and G',^ while flows corresponding 
to larger disorder (curves b, d, h) are affected only by the presence of FP G'(j. 

Fig. 121 presents the flow parameter dependencies of z^e for the curves in Fig. [T] The 
intrinsic feature of all curves is a non-monotonic behaviour of the effective dynamical 
critical exponent with approach to an asymptotic value. The experimental investigations 
are performed mainly in the non-asymptotic region. As it follows from Fig. |2l one 
can observe values of the dynamical exponent z that exceed or are smaller than an 
asymptotic one. Our asymptotic value z = 2.14 is somewhat smaller than the central 
value of the experimental result z = 2.18 ± 0.10 [SZj obtained for the dynamics of 
diluted Ising systems, however this experimental outcome is in agreement with our non- 
asymptotic observations. 




Figure 1. Static (left picture) and dynamic (right picture) flows of disordered model 
C for n = 1 {n < ric)- Filled circles mean unstable FPs with 7* = 0, blank circles 
indicate projections of unstable FPs with nonzero 7*, while the filled square denotes 
the stable FP. 

It is interesting to look at contributions of different origin to the effective dynamical 
critical exponent z^e- They are shown in Fig. Elfor flows a and h of FigQ correspondingly. 
The contributions consist of (i) the terms already present in model A (dashed curve in 
Fig. EJ, (ii) terms present in pure model C only (short dashed curve) and (iii) terms 
present in the diluted model C only (dashed-short dashed curve) [3]. The interplay of 
the above contributions gives the full effective exponent z^a (solid line) and may lead 
to an almost asymptotic value of the exponent although the parameters are far away 
from asymptotics. This is an important point, since the appearance of an asymptotic 
value in one physical quantity does not mean that other quantities have also reached 
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Figure 2. Dependencies of the effective dynamical exponent on the flow parameter 
for OP dimension n = 1 and n = 3. Curves are obtained on the basis of flows of Fig.Q 
and Fig. |3| 

the asymptotics. This is due to the different dependence of physical quantities on the 
model parameters. 




iog,„' iog,„' 

Figure 3. Different contributions to the effective dynamical critical exponent 
corresponding to the flows a and h of Fig. ^ The solid line represents the complete 
Zeff, for the other lines see the text. 

The second important case we consider here is the Heisenberg systems n = 3. 
In asymptotics their critical behaviour is characterized by the FP P. The static flows 
and projections of the dynamic flows in this case are shown in Fig. |3 There exists 
a considerable difference to the one-loop picture. In the one-loop approximation the 
stable FP P is marginal: for this FP every value of p* between and 1 is allowed. 
Therefore starting from different initial values flows approach different p* . In contrast, 
Fig. 13 shows that within two- loop approximation all flows remain near the initial values 
of p. Only when the flows reach the region m = m*, A = 0, p drops down to its FP value 
p* = 0. This may be attributed to the marginality in one-loop order. The behaviour of 
^eflf corresponding to flow a of Fig. |3] is shown in right picture of Fig. El We present the 
dependencies of the effective dynamical critical exponent on the flow parameter only for 
one set of values 7(^0) and p(^o) (but for two different ratios m(£o)/A(£o) ) because for 
other values 7(^0) and p{Iq) we observe similar behaviour: the peak increases for larger 
7(£o) and/or p(^o)- The experimental investigations give for a disordered Heisenberg 
magnet a value of the critical exponent z = 2.3 ± 0.1 [SH], that is larger than the value 



Model C critical dynamics of disordered magnets 



15 




Figure 4. Static (left picture) and dynamic (right picture) flows of disordered model 

C for n = 3 (n > Uc). Filled circles mean unstable FPs with 7* = 0, blank circles 
indicate projections of unstable FPs with nonzero 7*, while the filled square denotes 
the stable PP. 



5. Conclusions 

We have studied the critical dynamics of model C in the presence of disorder. 

Already one-loop order gives a qualitatively correct picture for suitable numbers 
of OP components n. In the asymptotics the conserved density is decoupled from the 
OP as it was expected. However the non-asymptotic critical behaviour of model C 
is strongly influenced by the presence of the static coupling between the OP and the 
conserved density as well as by the disorder. When approaching the asymptotic region 
the dynamical critical exponent usually shows a maximum starting from the van Hove 
value z = 2 in the background. For the case n < {uc = 4) the effective exponent 
might go after the maximum again through a minimum after it reaches its asymptotic 
value. It might be also the case that one observes the asymptotic value of the dynamical 
exponent although the parameters of the system are far from their asymptotic values. 
For systems with OP components n > ric only a maximum is observed. 

These general observations remain qualitatively the same in two-loop approximation 
and are modified only quantitatively. In particular, the value of the marginal dimension 
ric lies between 1 and 2 at d = 3. 

The existence of numerous FPs leads to a different cross-over behaviour of RG flows. 
Therefore the dynamical critical exponent can assume different values approaching 
asymptotics. For the case n < n^. z^s can cither exceed the asymptotic value showing a 
peak or reaches its asymptotic value mechanically. For a system with order parameter 
dimension n> z^s always shows a non-monotonous behaviour. 
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Appendix A. Dynamical functionals and perturbation expansion 

The model defined by expressions (P)-® within Bausch-Janssen- Wagner formulation 
[21] turns out to be described by an unrenormalized Lagrangian: 

d'^xdtl -f ^ (po,i(pQ,i+ ^ (po,i { — +f (r-V^) J v?o,i + A^moV^mo 

i=l i=l ^ ^ 

+faQ - a^AmV^ j mo + k ^ (^o.iV'o.i ^ V'ojV^oj 



+tV{x) J2 V?o,i(a;, t)ipo,i{x, t) + tj^mo ^ (po,npo,i 

i i 

-^Am7mmo VVo,iV?o,j. I (A.l) 

with new auxiliary response fields (pi{x,t). There are two ways to average over the 
random potential of disorder V{x) for dynamics. The first way is the same as in 
statics and consists in using the replica trick |,23j, where N replicas of the system 
are introduced in order to facilitate configurational averaging of the corresponding 
generating functional. Finally the limit — > has to be taken. 

However it was established j2S] that renormalization of the replicated Lagrangian 
leads to the same results as the renormalization of an Lagrangian obtained avoiding the 
replica trick but taking the mean of the Lagrangian ()A.1|) with respect to the random 
potential with distribution (jHJ. The Lagrangian obtained in this way can be written in 
the form 

£ = £0 + An* , (A.2) 

where the Gaussian part Cq is given by Q and Cint by (fTUj) . We perform the calculations 
on the basis of the Lagrangian defined by ()A.2|) using the Feynman graph technique. 
Response propagators for OP G{k,u) and secondary density H{k,u!) are equal to 

G{k,uj) = l/{-iuj + t{r + k'^)) and H{k,uj) = 1/ {-ico + UmX^k^) , (A.3) 

while the correlation propagators C{k,uj) and D{k,uj) are equal to 

= 2f/|-icj + f (r + P)|2 and L>(A;, cu) = 2A„A;V| -icu + a^A^fc^p .(A.4) 

We proceed within two-loop approximation. In order to obtain the two-point vertex 
function f V (f, A, 7^, f , A^, a^, k, uj) = f <^^(f, u, A, 7^, f , A^, a^, k, uj)5ij one needs 
to calculate diagrams of corresponding order. The result of the calculations can be 
expressed in the form: 

f ^^(e, k, u) = -tcuh^^i^, k, cu) + k)f. (A.5) 
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Here we introduce the correlation length ^(r, -u, A), which is defined by 



d In r^* 



17 



(A.6) 



fc2=0 



The function V^p^ is the static two- loop vertex function of a disordered magnet. The 
structure ()A.5|) of the dynamic vertex function of pure model C was obtained in |28j . 
We can express the two-loop dynamical function il(p^p in the form: 

where the one loop contribution has the structure 

1 



nJ^(e,fc,u;) = 4Af / 

Jk' 



-2cj + r(e-2 + A;'2))(e"2+A;'2) 
while the two-loop contribution is of the form: 



+ ^TIc{tk,uj), 



(A.7) 
(A.8) 
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4--^ruAW^';'''\^,k,u;)+16tA^W^';'''\^,^ 



4f (ly^J^^) i^,k,u;) + Wff'^ i^,k,uj) + 2W^( J^^) i^,k,uj)) , ( A.9) 

with the rescaled coupling 7 = 7m 

Expressions for the integral Ic and for ly'^"^^ C*'-^^-' and S of pure model C are given 
in Appendix A.l in Ref. ^UIj the contributions W^^^^'^ are the following: 




1 



'k' Jk" (^"' + + A;"2)(e-2 + (A; + A;' + A; 

1 



-icu + r(^-2 + (A; + A;' + A;")2)) 




r 



fc' Jk" + A;"2)(-icu + r(e-2 + A;'2))2(-zcj + 1(^-2 + k"^)) 



(A.IO) 
(A.ll) 



p|,(Ci?4) 




1 



k' Jk" + k'^){~tuj + r(e-2 + (A; + A;' + k"y)) 
1/1 f 



(^-2 + A;"2) \ (^"2 + + ^/ + p)2) 

£2 

-icu + f (e-2 + A;'2))(-zcj + f (^-2 + A;"2)) J 



-tcu + r(^-2 + k 



Il2\ 



(A.12) 



W, 



{CD5) 



r2 



fe'A"(^"'+^'")(-^^+r(e-'+A;'2))2 

1 

[-iuj+t{^-^+k"^)+x{k'+k"y) ' 



(A.13) 
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Ik' Jk" {i-^+k"^){-iu+v{i~^+k'^)+\{k+k'YY 
1 



'-iuj+v{i-^+k"^)+\{k+k'y) 




1 



'fc' Jk" + {k + k' + k"f){-iuj + r(e-2 + A;'2)) 
1 



f f2 



+ 



(A. 14) 



where we use rescaling A = amXm- 



(A.15) 



Appendix B. Renormalizing factors 

We use the minimal subtraction RG scheme [22] for renormahzation. In the definition 
of the renormahzing factors we follow Ref. [TU] . 

For renormahzation of the OP (p, fourth-order couplings m, A and correlation 
functions with insertion we introduce the renormalizing factors Z^p, Zu-, Z/^ and 
Z^p2 respectively via the relations: 

^o = ^yV, u = ii'Z~^Z^uA-,\ k = ii''^Z-^Z^^A-,\ l|y,o|2 = z^,l|y,|2^ (B.l) 

where is the scale and e = A — d. 

The renormahzation of r is performed via the relation 

r = Z~'Zrr. (B.2) 

The renormalizing factor Z^ is connected to Z^p2 by the relation Z^p2 = Z^^Zr therefore 
within minimal subtraction scheme one does not need to consider renormahzation for 
correlation functions containing ip'^ insertions explicitly. However a correlation function 



,2. 



containing two insertions < > needs additive renormahzation A,^ 

The renormahzation of dynamic quantities is introduced similar to statics. 
Renormalizing factors for the dynamic field ip, Z^, and kinetic coefficient F, Zr, are 
introduced via 

^,=zf^, t=ZrT. (B.3) 

The factor Z-p in the last equation contains a static contribution Z^^,, which can be 
separated: 

Zr = Z'J'Z^'^'zlf\ (B.4) 
Since in the dynamic model mode coupling terms are absent Z^^ = 1. Therefore 

Zr = Z]I^Z^'I\ (B.5) 
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In model C one needs to introduce additional renormalizing factors for the secondary 
density m and its coupling parameter 7. They are renormalized similar to and u in 
Eq. dEH]): 

ali^m, = Z^m, a-^^'^m = ^i''^ Z.^Af . (B.6) 

There are relations connecting the static Z-factors of model C to the static Z-factors 
of the Landau-Ginzburg- Wilson model by integrating out the secondary field m in the 
Hamiltonian (jSj). Thus the renormalizing factor of 7 is determined by: 

Z^ = Z^Z^Z^2. (B.7) 

Therefore Eq. ()B.6|1 can be rewritten as 

«;n'/'7m = /i^/' V^-7^d'^'- (B.8) 

The additive renormalization A^2 of the specific heat of the Landau-Ginzburg- Wilson 
model determines the renormalizing factor Z^ via the relation 



Z-2(m,A,7) = 1 + 7M^2(m). 



(B.9) 



Since the secondary density is conserved, no new renormalizing factor is needed for 
the dynamic auxiliary density m. It renormalizes by: 

Z~^rh. 



The kinetic coefficient A renormalizes as 
O'm^m = Z\X. 

Similar to Z^ its static contribution can be separated: 



72 7 



(d) 



Since no mode coupling terms are present in model C, zi^"^ = 1, therefore 



(B.IO) 
(B.ll) 
(B.12) 
(B.13) 



Renormalizing T^^ we obtain the two-loop renormalizing factor Z^ in the form: 

„2 „„ 1 . / / I ^\\ n + 2 



w 



A 7 
1-8- -2I 

£ e 1 + w 



1 

+ ^ 



w 



7 



1 + w \ w + 1 



+4A72 

+7 



w 



1+w 
2 w (n 



+2OA72 



w 



n 

- - 1 
2 



n+2 



w 



u 7 



+ 96A2-4^-^mA 



+ 



2e 



l+w\2 1+^ 2 3~ 1+w ^ l+2w 



7 



n+2 

2 w 



u 1-3 In- 



\+w 



n+2 , 
+3— ^mA 



-44A2-I2A72 



w 



\+w 



n+2 u 



12 e 



- In--- + 



3 12 



2 ""^ 



e 1 + w 



w 



w ^ w 3 , ,1 

— m m 1 + u; 

2 1+w 2 ^ ' 21+w 



\nw 



(B.14) 
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